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The  executive  summary  encapsulates  the  work  of  the  AFOSR/NASA  National  Science  Center  for 
Research  in  Hypersonic  Laminar-Turbulent  Transition  hereafter  referred  to  as  NHSC.  It  includes 
a  brief  description  of  accomplishments  in  the  order:  (1)  Executive  Summary;  (2)  Experimental 
Eacilities  and  Results;  (3)  Theoretical  and  Computational  Modeling;  (4)  Students  and  Publications. 
This  is  followed  by  ten  volumes  of  details  contained  in  dissertations. 

1.  EXECUTIVE  SUMMARY 

1.1  Overview 

The  NHSC  was  originally  established  for  5  years  beginning  1  April  2009  at  a  funding  level  of  $2 
million  per  year  with  AEOSR  paying  the  first  half  and  NASA  paying  the  second  half  Dr.  William 
Saric  was  the  Center  Director.  Due  to  sequestration  in  the  spring  of  2013,  NASA  was  unable  to 
honor  its  last  increment  of  funding  so  that  it  was  effectively  a  four  and  one-half  year  program. 
Since  the  Texas  A&M  Board  of  Regents  approved  the  “Center”,  it  was  decided  to  keep  this 
designation  for  future  funding  and  to  have  Drs.  Rodney  Bowersox  and  Helen  Reed  as  co-Directors. 

The  NHSC  was  made  up  of  1 1  Professors  and  2  consultants.  The  participating  universities  were 
California  Institute  of  Technology  (Caltech),  Texas  A&M  University  (TAMU),  University  of 
California  -  Eos  Angeles  (UCLA)  and  University  of  Arizona  (UA).  TAMU  was  the  lead  and 
subcontracted  to  the  other  universities  and  consultants.  The  participants  were  chosen  from  the 
leaders  of  the  high-speed  transition  community  who  covered  all  of  the  relevant  areas  of  analysis, 
computation,  experiment,  facilities,  and  diagnostics.  The  group  had  one  dedicated  Technical 
Exchange  Meeting  (TEM)  each  year  whose  venue  rotated  over  the  four  universities.  They  also 
conducted  mini-TEMs  twice  a  year  at  the  two  relevant  AIAA  meetings  (Jan/June).  These  were 
open  to  the  public. 

In  this  summary,  the  format  [**]  is  used  to  reference  the  relevant  publication  that  was  produced  as 
part  of  the  center.  Eollowing  the  Executive  Summary,  the  detailed  reports  from  each  co-PI  are 
contained  in  individual  volumes. 

1.2  Highlights 

1.2.1  Three  major  hypersonic  facilities  made  operational  at  TAMU  and  a  fourth,  T5  at  Caltech, 
was  improved  so  that  it  was  able  to  measure  Mack  modes. 

1 .2.2  Three  new  diagnostic  tools  developed  (EST,  ELDI,  VENOM). 

1.2.3  Nonlinearities  measured  with  EST.  Then,  Klebanoff-type  streaks  were  observed  with  IR 
Thermography.  This  complements  the  computational  predictions  of  Easel  (3.4)  and  Reed 
(3.1). 

1 .2.4  Eirst  measured  transient  growth  streaks  in  a  hypersonic  boundary  layer. 

1.2.5  Developed  in-house  LST/PSE/NPSE  with  V&V  for  comparisons  with  experiments  and 
with  DNS  (nonlinear  breakdown). 

1 .2.6  Explained  the  pressure-based  mechanisms  that  lead  to  stabilization  in  hypersonic  flows. 

1.2.7  Showed  that  pressure-field  evolution  exhibits  a  three-stage  transition  process  which  was 
also  shown  in  the  DNS. 

1 .2.8  Clarified  the  terminology  for  the  second-mode. 

1 .2.9  Demonstrated  that  discrete-mode  branching  can  be  resolved  with  PSE. 
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1 .2. 10  The  second-mode  stabilization  with  porous  coatings  has  been  put  on  firm  ground  with  the 
combination  of  theory,  DNS,  and  experiment  done  by  Fedorov,  Homung,  Fasel,  and  Zhong 
as  part  of  the  NHSC. 

1.2.11  As  another  example  of  multiple  Pis  interacting  across  the  center,  useful  models  for 
including  chemistry  have  been  advanced  by  Tumin,  Girimaji,  Fasel,  Reed,  and  Zhong  as 
part  of  the  NHSC. 

1.2.12  The  spatial  DNS  calculations  showed  that  nonlinear  wave  interactions  (Klebanoff-type 
breakdown  with  streaks),  the  Gaster  wave  packet,  and  the  turbulent  spot  all  had  remarkable 
similarities  to  what  was  observed  in  incompressible  boundary  layers.  Thus,  it  appears  that 
we  already  can  classify  the  fundamental  mechanisms  for  breakdown  in  hypersonic 
boundary  layers  and  we  need  to  concentrate  on  the  details. 

1.2.13  41  PhD  students  were  supported  by  NHSC.  47  journal  articles  and  130  conference 
proceedings  were  written  by  NHSC. 

1.3  Recommendations 

The  Orr-Sommerfeld  Equation  (OSE)  is  1 1 1  years  old  and  progress  in  boundary-layer  stability  and 
transition  has  been  slow  -  much  of  it  occurring  over  the  past  50  years.  A  review  of  the  history 
shows  that  major  breakthroughs  have  been  made  on  a  scale  of  10  years.  This  doesn’t  mean  that 
once  every  ten  years  someone  has  a  good  idea.  Rather,  given  an  area  of  stability/transition  research 
such  as  receptivity,  nonlinear  interactions,  3-D  boundary  layers,  transient  growth,  etc.,  it  takes  ten 
years  of  continuous  hard  work  on  a  transition  problem  in  order  to  make  a  meaningful  contribution. 

The  team  structure  of  the  Center  really  came  together  after  the  infrastructure  was  established. 
Productivity  and  cooperation  increased  each  year  and  important  results  were  being  generated.  This 
was  a  worthwhile  activity  for  everyone  associated  with  the  Center.  If  transition  continues  to  be  an 
important  issue  for  the  Air  Eorce,  a  successful  working  model  has  been  established.  It  took  a  heroic 
effort  to  support  the  multiple  PI  NHSC  and  it  was  difficult  to  sustain  the  funding.  Perhaps  small 
computational/experimental  teams  can  be  supported  out  of  the  core  program. 

Before  things  get  too  far  along,  we  should  like  to  suggest  an  AFOSR/NASA-sponsored 
international  workshop  on  hypersonic  transition  that  can  build  on  the  very  successful  NHSC 
supported  workshop  held  in  Sedona  2012.  The  NHSC  can  then  review  all  of  their  accomplishments 
for  the  entire  community.  It  was  thought  that  the  NCTAM  held  this  year  would  do  this  but  the 
conflict  with  the  AIAA  summer  meeting  prevented  this  from  happening. 

New  ideas  are  necessary.  Approaching  receptivity  and  stability  as  an  initial-value  problem  is  the 
only  way  to  get  away  from  eigenvalue  calculations.  Solutions  of  OSE  only  provide  correlations. 
The  federal  agencies  must  be  willing  to  support  speculative  diagnostic  development.  Amplitude 
growth  within  the  boundary  layer  is  O(IO^)  and  we  need  to  quantify  the  freestream  where 
disturbances  in  the  unstable  pass  band  may  be  0(10'^).  We  have  to  continue  to  develop  the  ideas 
of  Eedorov  [60]  that  look  at  not  only  thermal  fluctuations  as  a  least  upper  bound  on  freestream 
disturbances  but  freestream  vorticity  and  surface  roughness  as  well.  These  latter  two  have  a 
profound  effect  on  3-D  boundary  layers. 

Transition  is  a  difficult  problem  and  support  must  be  continuous.  The  students  that  are  graduated 
must  be  encouraged  and  supported  to  keep  working  in  the  area. 
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2.  EXPERIMENTAL  FACILITIES,  DIAGNOSTICS,  AND  RESULTS 

The  objective  of  the  foundational  experiments  was  to  perform  fundamental  experiments 
coordinated  with  theoretical  and  numerical  requirements.  The  outcome  of  this  work  was  to  provide 
a  comprehensive  experimental  map  of  the  hypersonic  transition  parameter  space  and  establish  a 
database  of  mechanisms  and  growth  rates  in  quiet  flow.  Considerable  effort  was  also  placed  on 
facility  and  instrumentation  development.  The  experimental  accomplishments  are  listed  below: 

(1)  Facility  and  Instrumentation  Advancement 

a.  Re-establishment  and  calibration  of  the  Mach  6  Quiet  Tunnel  at  TAMU 

b.  Incorporation  of  novel  optical  diagnostics  in  the  calibration  of  the  Caltech  T5  tunnel 

c.  Calibration  of  the  TAMU  ACE  Hypersonic  Tunnel 

d.  Further  development  of  the  laser-based  VENOM  diagnostic 

e.  Design  and  construction  of  the  TAMU  Mach  5  shock  tunnel 

f  Development  of  the  focused  schlieren  technique  (FST)  and  implementation  for 
nonlinear  measurements  and  angle-of-attack  alignment  at  TAMU. 

g.  Development  of  the  Focused  Easer  Differential  Interferometry  (FEDI)  at  Caltech. 

h.  Design  and  fabrication  of  a  blunt-body  model 

(2)  Hypersonic  Stability  and  Transition  Measurements 

a.  Second  mode  instability  on  the  Langley  93-10  cone  in  the  M6QT 

b.  Second  model  instability  under  high-enthalpy  conditions  in  T5 

c.  Transient  growth  instability  on  a  5-degree  cone  in  the  M6QT 

d.  Trip  instability  and  transition  in  the  ACE  Tunnel 

e.  Crosshatch  roughness  and  pressure  gradient  effects  in  the  TAMU  High  Re  tunnel 

f  Characterize  the  freestream  and  boundary  layer  in  T5  using  FEDI. 

g.  Apply  the  FST  with  bi-coherence  to  identify  the  higher  harmonics  in  the  transitional 
region  as  nonlinear  interactions  of  the  fundamental. 

Collectively,  the  experimental  efforts  resulted  in  new  experimental  capabilities  and  carefully 
acquired  databases  for  theoretical  model  development  and  validation.  It  is  expected  that  these 
diagnostics  will  be  adapted  for  use  in  other  hypersonic  facilities.  A  brief  synopsis  of  the  role  of 
each  PI  is  given  below.  Detailed  descriptions  are  given  in  the  attached  volumes. 

2,1  Quiet  Tunnel  Experiments  (William  Saric):  Experiments  relevant  to  flight  require  ground- 
test  facilities  with  very  low  disturbance  levels.  To  enable  this  work,  the  Mach  6  Quiet  Wind  Tunnel 
(M6QT),  previously  of  NASA  Langley  Research  Center,  was  reestablished  within  a  new  pressure- 
vacuum  blowdown  infrastructure  at  Texas  A&M.  A  40-second  run-time  at  constant  conditions 
enables  detailed  measurements  for  comparison  with  computation.  The  freestream  environment 
was  extensively  characterized,  with  a  large  region  of  low-disturbance  flow  (0.05%)  found  to  be 
reliably  present  for  unit  Reynolds  numbers  Re'  up  to  1 1  x  10^  /m  [76] 

The  first  set  of  quiet  flow  experiments  were  performed  by  the  Saric  group  on  a  5“  half-angle  flared 
cone  model  at  Re'  =  10x10^  /m  and  zero  angle  of  attack.  For  the  study  of  the  second-mode 
instability,  well-resolved  boundary-layer  profiles  of  mean  and  fluctuating  mass  flux  were  acquired 
at  several  axial  locations  using  hot-wire  probes  with  a  bandwidth  of  330  kHz.  The  second  mode 
instability  was  observed  to  undergo  significant  growth  between  250  and  310  kHz.  Mode  shapes  of 
the  disturbance  agree  well  with  those  predicted  from  linear  parabolized  stability  equation  (LPSE) 
computations.  A  17%  (40  kHz)  disagreement  is  observed  in  the  frequency  for  most-amplified 
growth  between  experiment  and  LPSE.  Possible  sources  of  the  disagreement  are  discussed  in  the 
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detailed  volumes.  The  effect  of  small  misalignments  of  the  model  is  quantified  experimentally.  A 
focused  schlieren  deflectometer  with  high  bandwidth  (1  MHz)  and  high  signal-to-noise  ratio  was 
employed  to  complement  the  hot-wire  work.  The  second-mode  fundamental  at  250  kHz  was 
observed,  as  well  as  additional  harmonic  content  not  discernible  in  the  hot-wire  measurements  at 
two-  and  three-times  the  fundamental.  A  bi-spectral  analysis  showed  that  after  sufficient 
amplification  of  the  second  mode,  several  nonlinear  mechanisms  become  significant,  including 
ones  involving  the  third  harmonic,  which  have  not  hitherto  been  reported  in  the  literature  [13,  76, 
77,  78,  79,  80]. 

The  are  some  nascent  projects  in  the  M6QT  involving:  (1)  the  7°  yawed  straight  cone  for  crossflow 
to  go  along  with  Reed’s  computations  [127,  132];  (2)  the  5°  flared  cone  examine  Klebanoff-type 
breakdown  in  transition  as  demonstrated  by  Fasel  and  Reed;  (3)  the  cooled  cone  from  Purdue  to 
enhance  transition.  This  work  is  moving  at  a  snail’s  pace  since  sequestration. 

The  details  of  this  work  are  contained  in  the  Volume  V  supplement. 

2.2  Diagnostics  and  Roughness  Experiments  (Rodney  Bowersox):  In  addition  to  collaborating 
with  the  Saric  group  on  the  M6QT  installation  and  the  North  group  on  the  VENOM  diagnostic, 
the  Bowersox  group  developed  and  calibrated  the  Actively  Controlled  Expansion  Hypersonic 
Wind  Tunnel  (ACE),  which  is  a  0.23  x  0.36  m,  variable  Mach  (5  -  8)  facility.  This  tunnel  was  used 
for  late  stability  and  transitional  measurements.  The  freestream  mean  flow  uniformity  was  less 
than  0.5%.  The  freestream  disturbance  levels  were  found  to  vary  from  0.15%  to  1.5%  depending 
Reynolds  number.  Also  a  new  Mach  5  shock  tunnel  (0.13  m  x  0.13  m  nozzle  exit)  was  developed 
to  aid  in  porting  the  VENOM  diagnostic  to  high-enthalpy  impulse  tunnels.  Measurements  were 
also  made  in  the  supersonic  high-Reynolds-number  tunnel  operating  at  Mach  5  [87,  141,  143,  145, 
146,  158]. 

NASA  diamond  (“pizza  box”)  trips  were  used  within  a  low-Reynolds-number  boundary  layer  (Reo 
=  3700),  which  allows  for  DNS  computations.  Hot-wire  and  Pitot  pressure  measurements  were 
acquired.  Moving  downstream  of  the  trips,  the  peak  RMS  disturbance  value  grew  in  amplitude 
until  saturation.  After  this  point,  the  peak  decreased  in  amplitude,  spread  throughout  the  boundary 
layer,  and  ultimately  led  to  transition.  Detailed  measurements  were  then  acquired  in  the  late 
transitional/turbulent  boundary  layer.  The  power  spectra  followed  the  traditional  power  law,  which 
verified  that  the  length  scales  followed  the  behavior  described  by  Kolmogorov.  The  profile  was 
quantified  across  the  boundary  layer  into  the  sublayer  [17,  40,  156,  157]. 

The  effects  of  crosshatch  roughness  {k/S  ~  0.07)  and  pressure  gradient  {Ip  =  -0.08  and  -0.49)  on  a 
Mach  5  high  Reynolds  number  {Ree  =  40,000)  turbulent  boundary  layer  were  quantified  using 
high-resolution  particle  image  velocimetry  (PIV).  The  distortion  and  reorientation  of  the  large- 
scale  coherent  motions  is  quantified  through  the  determination  of  the  integral  length  scale  and  local 
structure  angle  from  two-point  correlations.  Detection  of  individual  hairpin  vortices  through  the 
swirling  strength  criterion  Id  allows  the  population  distribution  of  the  turbulent  eddies  to  be 
examined,  along  with  the  conditionally  averaged  hairpin  structure.  The  reduced  Reynolds  stresses 
observed  in  the  favorable  pressure  gradients  is  partially  due  to  the  attenuation  of  the  local  fiowfield 
around  the  near-wall  hairpin  structures,  mitigating  the  mechanism  for  “producing”  turbulence.  The 
rotational  rate  of  the  hairpin  vortices,  measured  through  the  mean  prograde  swirling  strength,  was 
reduced  for  the  favorable  pressure  gradient  models  [2,  87,  122,  123,  124,  125]. 

The  details  of  this  work  are  contained  in  the  Volume  I  supplement. 
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2.3  Blunt-Body  Experiments  (Eli  Reshotko  with  Bowersox,  Hornung,  and  Reed):  A  blunt- 
body  model,  in  the  shape  of  the  MSL,  was  fabrieated  for  determining  the  effeet  of  roughness  at 
angle  of  attaek.  Bloekage  tests  were  suceessfully  completed  in  ACE.  This  work  was  scheduled  for 
the  last  half-year  of  the  center  but  as  mentioned  earlier,  the  project  was  terminated.  [93,  110,  258, 
260,  263,  264] 

2.4  Diagnostics  (Simon  North):  The  North  group  worked  to  advance  the  Vibrationally  Excited 
Nitric  Oxide  Monitory  (VENOM)  diagnostics,  which  combines  two-line  Planar  Easer  Induced 
Eluorescence  methods  and  two-component  Molecular  Tagging  Velocimetry.  The  technique 
monitors  the  nascent  NO(v"=l)  arising  from  photo-dissociation  of  trace  amounts  of  NO2  as  a 
molecular  tracer.  The  VENOM  technique  is  expected  to  be  not  only  applicable  to  cold  high-speed 
flows,  which  is  the  focus  of  the  present  work,  but  also  to  combustion  and  other  reactive  or  high- 
enthalpy  flow  fields.  During  this  program,  the  spatial  resolution  was  refined  to  0.3  mm,  with  a 
velocity  and  temperature  uncertainty  of  1%  and  5%,  respectively.  In  addition,  a  new  hypersonic 
flow  apparatus  was  developed  for  advanced  laser  diagnostics  development.  This  apparatus  is 
characterized  by  its  pulsed  operation  mode  that  translates  into  a  significant  reduction  in  mass  flow 
rates  and  can  be  operated  for  long  periods  at  Mach  numbers  ranging  from  2.8  to  6.2.  The  flow 
conditions  during  the  uniform  flow  time  interval  of  each  pulse  vary  by  less  than  1%,  generating  a 
flow  of  sufficient  quality  for  quantitative  measurements  [15,  16,  32,  33,  34,  83]. 

The  details  of  this  work  are  contained  in  the  Volume  VI  supplement. 

2.5  Transient  Growth  Experiments  (Edward  White):  The  White  group  performed  the  second 
major  experiment  in  the  M6QT.  This  study  focused  on  quantifying  the  effects  of  surface  roughness 
on  boundary-layer  disturbance  growth  and  laminar-to-turbulent  transition.  The  transient  growth 
mechanism  that  produces  algebraic  growth  of  streamwise  streaks  via  decaying  streamwise  vortices 
has  not  previously  been  deliberately  observed  in  hypersonic  flow.  The  measurements  were 
performed  on  in  the  boundary  layer  on  a  5°  half-angle  smooth  cone  fitted  with  a  slightly  blunted 
nosetip  and  a  ring  of  18  periodically-spaced  cube-like  discrete  roughness  elements  1-mm  tall  by 
1.78-mm  wide  by  1.78-mm  long.  The  roughness  element  height  was  approximately  equal  to  the 
boundary-layer  thickness,  yet  no  transition  to  turbulence  is  observed  for  freestream  unit  Reynolds 
numbers  between  7.5  x  10^  /m  and  9.8  x  /m.  Pitot  measurements  revealed  azimuthally- 
alternating  high-  and  low-speed  streaks  growing  downstream  of  the  roughness.  Earge  unsteadiness 
was  measured  in  the  roughness  wake  but  decays  downstream.  The  streamwise  evolution  of  the 
steady  and  unsteady  disturbance  energy  is  consistent  with  low-speed  observations  of  transient 
growth  in  the  mid-wake  region  behind  periodically-spaced  cylindrical  roughness  elements.  This 
experiment  contains  the  first  quantitative  measurements  of  roughness-induced  transient  growth  in 
a  high-speed  boundary  layer  [111,  147]. 

The  details  of  this  work  are  contained  in  the  Volume  IX  supplement. 

2.6  T5  Experiments  at  Caltech  (Joseph  Shepherd  and  Hans  Hornung):  The  Caltech  team 
focused  on  incorporating  novel  optical  techniques  to  quantify  the  freestream  disturbances  in  T5 
and  second  mode  instability  growth  on  a  slender-body  under  hypervelocity  conditions,  which  is  a 
previously  unexplored  regime  where  thermochemical  effects  are  important  [14].  Narrowband 
disturbances  (0. 5-3.0  MHz)  are  measured  in  boundary  layers  with  edge  velocities  of  up  to  5  km/s 
at  two  points  along  the  generator  of  a  5  degree  half  angle  cone.  The  freestream  disturbances  were 
found  to  be  order  0.1%  in  the  second-mode  passband  [12,  28].  Experimental  amplification  factor 
spectra  were  acquired.  Einear  stability  and  PSE  analysis  was  performed,  with  fair  prediction  of  the 
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frequency  content  of  the  disturbances;  however,  the  analysis  over-predicts  the  amplification  of 
disturbances.  The  results  of  this  work  have  two  key  implications:  1)  the  acoustic  instability  is 
present  and  may  be  studied  in  a  large-scale  hypervelocity  refiected-shock  tunnel,  and  2)  these  data 
provide  a  new  basis  on  which  the  instability  can  be  studied  [45,  81,  82,  116,  120]. 

The  development  of  the  focused  laser  differential  interferometry  (FLDI)  technique  was  certainly 
one  of  the  highlights  of  the  center  [27,  117,  118,  121].  This  permitted  high  bandwidth 
measurements  in  both  the  freestream  and  the  boundary  layer  of  T5  with  the  result  that  there  is  very 
little  spectral  content  in  the  freestream  disturbances  in  the  range  of  1  -  2  MHz.  This  is  the  second¬ 
mode  (Mack  mode)  range  for  high  enthalpy  (high  speed)  boundary  layers.  The  implications  are 
huge  in  that  this  (amongst  other  things  of  course)  probably  accounts  for  the  success  of  the  porous 
coatings  laminar  flow  control  [3,  4,  8]. Laminar  flow  control  using  C02  injection  was  also  studied. 
[89,  90,  91,  92,  105,  164].  It  is  fair  to  say  that  these  are  promising  areas  of  further  study. 

The  details  of  this  work  are  contained  in  the  Volume  VII  supplement. 

3.  THEORETICAL  AND  COMPUTATIONAL  MODELLING 

Over  the  past  decades  of  studying  stability  and  transition,  it  has  become  apparent  that  it  is  critically 
important  for  computations  and  experiments  to  work  very  closely  together  on  the  same  geometries 
and  operating  conditions.  Advances  in  basic  mechanisms  and  prediction  methods  have  come  from 
working  together.  Transition  is  highly  sensitive  to  operating  conditions,  especially  in  the 
hypersonic  range.  Computations  provide  validation  of  experiments  and  vice  versa  [30,  134,  135, 
136,  137]. 

3.1  NPSE  Calculations  -  Verification  &  Validation  (Helen  Reed):  Within  the  NHSC,  the 
computational  group  of  Helen  Reed  has  been  collaborating  with  the  other  investigators  and 
providing  two-way  verification  and  validation  of  findings,  and  a  bridge  from  receptivity  to  the 
nonlinear  stages  of  transition.  The  objectives  were  to  extend  the  existing  analytical  framework 
beyond  Mack’s  linear  stability  formulation  to  include  relevant  3-D  hypersonic  flow  physics, 
including  identifying  the  relevant  instabilities  and  interactions,  carefully  validating  the  predictive 
tools  created,  and  disseminating  these  tools  once  validated.  This  group  has  been  modeling 
attachment-line  problems.  Mack-mode  instabilities,  and  crossflow  instabilities,  including  on 
yawed/unyawed  straight  and  flared  cones,  and  more  recently  on  elliptic  cones. 

Part  1:  Attachment-line  heating/cooling  in  a  compressible  flow  [31,  133]  was  studied  for  mainly 
laminar  flow  control  issues  with  the  idea  of  either  enhancing  or  ameliorating  one  (little  studied) 
fundamental  instability  on  swept  wings  and  elliptic  cones.  Only  modest  temperature  differences 
are  needed  to  affect  this  instability. 

Part  2:  JoKHeR:  NPSE  simulations  of  hypersonic  crossflow  instability  were  developed  [97,  98, 
126,  131].  It  was  immediately  apparent  that  an  in-house  code  was  required  for  prediction  of 
unstable  waves  with  inhomogeneous  boundary  conditions.  This  code  has  been  used  for  virtually 
all  of  the  TAMU  stability  calculations.  It  also  supports  the  on-going  experimental  work  in  the 
M6QT  on  nonlinear  breakdown  of  Mack  modes. 

Part  3:  Instabilities  on  a  T  hypersonic  yawed  straight  cone  are  concerned  principally  with 
crossflow  instabilities  and  this  work  supports  the  nascent  experimental  work  in  the  M6QT  [127, 
132]. 

Part  4:  Boundary-layer  instability  and  transition  on  a  flared  cone  [96]. 
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Part  4a:  Boundary-layer  instability  and  transition  on  a  flared  cone  in  a  Mach  6  quiet  wind  tunnel 
[79] 

Part  4b:  Stability  of  hypersonic  compression  cones  [126]. 

Part  4c:  Hypersonic  stability  analysis  of  a  flared  cone  [96]. 

Part  5:  Verification  and  validation  Issues  in  hypersonic  stability  and  transition  prediction  [134]. 
The  details  of  this  work  are  contained  in  the  Volume  IV  supplement. 

3.2  Transition  Modelling  (Sharath  Girimaji):  The  objectives  of  the  Girimaji  group’s  research 
was  to  incorporate  their  unique  research  in  scalar  transport  in  turbulent  flows  to  later  stages  of 
boundary  layer  transition  with  the  idea  of  bridging  backward  from  low-Reynolds-number 
turbulence  to  the  later  stages  of  transition. 

Part  la:  Instability  mechanisms  in  hypersonic  shear  flows.  In  incompressible  flows,  pressure  is  a 
Lagrange  multiplier  with  the  sole  purpose  of  imposing  the  divergence-free  constraint  on  the 
velocity  field.  In  high  speed  flow,  pressure  is  a  bona  fide  thermodynamic  variable  that  evolves 
according  to  a  wave  equation  derived  from  conservation  of  energy  and  state  equations.  This  wave 
behavior  of  pressure  introduces  new  flow  mechanisms  that  strongly  influence  flow  stability, 
transition  and  ultimately  turbulence.  The  objective  of  this  portion  of  work  is  to  examine  and 
explain  the  pressure-based  mechanisms  that  lead  to  flow  stabilization  in  hypersonic  flows  [36,  37, 
38,46]. 

Part  lb:  Three-stage  Instability  Mechanism:  Examination  of  the  linearized  perturbation  equations 
reveals  that  the  pressure-field  evolution  exhibits  three  distinct  regimes  in  a  hypersonic  shear  flow. 
In  the  first  regime,  the  timescale  of  shear  is  much  larger  than  that  of  pressure  evolution.  This  leads 
to  perturbation  velocity  field  evolution  purely  according  to  the  dictates  of  production  mechanism. 
The  second  regime  represents  the  state  in  which  the  shear  and  pressure  timescales  are  comparable. 
Under  these  conditions,  the  shear-normal  velocity  and  pressure  fields  interact  as  in  a  harmonic 
oscillator  and  yield  and  negative  production.  The  production  is  nearly  zero  when  integrated  over 
each  acoustic  cycle  causing  the  perturbations  to  be  stabilized.  In  the  third  regime,  the  acoustic 
timescale  is  smaller  than  that  of  shear  and  the  instabilities  evolve  in  a  manner  similar  to  that  in 
incompressible  turbulence  [1,  23]. 

Part  Ic:  Perturbation  Obliqueness  effects:  Velocity-pressure  dynamics  of  individual 
perturbation/fiuctuation  modes  is  investigated  using  direct  numerical  simulations  and  linear 
analysis  in  high  Mach  number  homogeneous  shear  flow.  For  a  given  perturbation  wave  mode,  the 
action  of  pressure  is  shown  to  depend  on  two  important  factors:  the  orientation  of  the  wave  vector 
with  respect  to  the  shear  plane  and  modal  Mach  number.  It  is  shown  that  the  streamwise 
perturbation  wave  modes  rapidly  develop  high  level  of  dilatation  but  are  self-limiting  due  to  the 
action  of  pressure.  On  the  other  hand,  spanwise  perturbation  wave  modes  develop  purely 
solenoidal  velocity  fluctuations  and  are  unaffected  by  pressure  or  Mach  number.  Oblique  modes 
(oblique  orientation  with  shear  plane)  combine  solenoidal  and  dilatational  characters  and  are 
shown  to  be  chiefly  responsible  for  stabilizing  effects  seen  in  engineering  flows.  Three  regimes  of 
obliqueness  of  different  stability  characteristics  are  identified.  The  investigation  also  isolates  the 
linear,  non-linear,  inertial  and  pressure  contributions  in  an  attempt  toward  a  comprehensive 
explanation  of  compressibility  effects  in  high  Mach  number  shear  flows.  The  effects  of  linear  and 
quadratic  velocity  profiles  have  been  examined  in  detail. 
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Part  Id:  Possible  instability  control  mechanism:  The  key  role  played  by  pressure  in  compressible 
flows  presents  new  avenues  for  the  formulation  of  novel  control  strategies  for  controlling 
instability  growth.  Such  control  strategies  can  have  foundations  in  the  energy  equation  and  exploit 
the  strong  coupling  between  momentum  and  energy  balance  equations  It  is  demonstrated  that 
'equi-partition’  of  kinetic  and  potential  energy'  observed  in  second  stage  of  growth  of  kinetic 
energy  observed  in  sheared  flow  can  be  exploited  to  bring  about  the  flow  control.  The  so-called 
potential  energy  of  turbulence  resides  in  the  fluctuating  dilatational  pressure  and  exchanges  energy 
with  dilatational  kinetic  energy  through  the  pressure  dilatation  term.  If  we  can  devise  control 
mechanisms  to  drain  energy  from  potential  energy  through  thermodynamic  or  acoustic  means,  then 
the  level  of  wall  normal  velocity  will  also  decrease  signiflcantly.  Since  the  production  of 
turbulence  in  shear  flows  is  highly  dependent  on  this  component  of  velocity  fluctuation,  the 
subsequent  growth  of  fluctuations  will  be  substantially  reduced. 

Part  le:  Breakdown  toward  turbulence:  It  is  demonstrated  that  breakdown  toward  turbulence  is 
signiflcantly  delayed  in  hypersonic  shear  flows  due  to  the  suppression  of  spanwise  vortical 
structures.  This  effect  is  brought  about  by  the  action  of  pressure  in  Regime-2  as  described  above. 
It  is  shown  that  non-linear  interactions  produce  both  streamwise  and  spanwise  vortical  structures 
at  higher  harmonics.  The  streamwise  vortical  structures  grow  unimpeded  by  compressibility 
effects  but  cannot  by  themselves  lead  to  breakdown  toward  turbulence.  It  is  shown  that  at  later 
times,  high-wavenumber  low-Mach  number  spanwise  vortical  structures  appear  due  to  non-linear 
effects.  The  combination  of  these  high-harmonic  spanwise  vortices  and  streamwise  vortices 
ultimately  leads  to  breakdown  toward  turbulence. 

Part  2:  Gas  Kinetic  Scheme  for  hypersonic  non-equilibrium  flows 

Part  2a:  Gas  kinetics-based  Numerical  scheme  for  hypersonic  non-equilibrium  flows:  The  recently 
developed  Gas  Kinetic  Method  (GKM)  for  fluid  flow  computations  is  enhanced  with  advanced 
reconstruction  (interpolation)  schemes  to  enable  direct  simulation  of  highly  compressible  turbulent 
fields.  Variants  of  Weighted  essentially  non-oscillatory  (WENO)  reconstruction  schemes  of 
different  orders  of  accuracy  are  implemented  and  examined  along  with  the  more  elementary  van 
Leer  method.  The  different  schemes  are  evaluated  for  their  accuracy,  efficiency  and  numerical 
stability.  The  computed  results  are  compared  against  the  rapid  distortion  theory  (RDT)  for  the  case 
of  compressible  shear  turbulence  and  'pressure-released'  Burgers  solution  at  high  enough  Mach 
numbers.  In  the  case  of  decaying  isotropic  turbulence,  the  efficacy  of  the  reconstruction  schemes 
is  evaluated  by  comparison  against  a  gold  standard'  high  resolution  simulation.  The  capabilities 
of  the  reconstruction  schemes  to  capture  linear,  non-linear,  pressure-released  and  viscous  flow 
physics  as  well  as  solenoidal  and  dilatational  features  of  the  flow  fields  are  established  in  isolation 
and  combination.  The  most  suitable  WENO  variant  for  integration  with  GKM  is  identified. 
Another  important  outcome  of  the  study  is  the  finding  that  temperature-interpolation  is  superior  to 
energy-interpolation  during  the  process  of  information  transfer  from  cell-center  to  cell-interface. 
Overall,  this  work  advances  the  applicability  of  kinetic  theory  based  GKM  to  a  wider  range  of  high 
Mach  number  flow  physics,  specifically  hypersonic  boundary  layers  [20,  23,  46]. 

Part  3:  Reduced  chemical  kinetics  modeling  for  hypersonic  flow  applications 

In  external  hypersonic  flows,  viscous  and  compressibility  effects  generate  very  high  temperatures 
leading  to  significant  chemical  reactions  among  air  constituents.  Therefore,  hypersonic  flow 
computations  require  coupled  calculations  of  flow  and  chemistry.  Accurate  and  efficient 
computations  of  air-chemistry  kinetics  are  of  much  importance  for  many  practical  applications  but 
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calculations  accounting  for  detailed  ehemical  kineties  can  be  prohibitively  expensive.  In  this  work, 
we  investigate  the  possibility  of  applying  ehemical  kinetics  reduction  schemes  for  hypersonic  air- 
chemistry.  We  consider  two  chemieal  kinetics  sets  appropriate  for  three  different  temperature 
ranges:  2500K  to  4500K;  4500K  to  9000K;  and  above  9000K.  By  demonstrating  the  existenee  of 
the  so-called  the  slow  manifold  in  each  of  the  ehemistry  sets,  we  show  that  judicious  chemical 
kinetics  reduction  leading  to  significant  computational  savings  is  possible  without  mueh  loss  in 
aeeuracy  [11,  88]. 

The  details  of  this  work  are  eontained  in  the  Volume  III  supplement. 

3,3  Foundational  Analysis  and  Control  (Anatoli  Tumin  and  Alexander  Fedorov):  As  part  of 
the  physics-based  approach,  it  was  recognized  early  on  that  a  fundamental  foundation  was 
necessary  for  the  understanding  of  the  “second-mode”  (or  Maek  mode)  and  as  well  as  the  role  of 
chemistry.  Analytical  and  computational  approaches  were  used  to  address  these  issues. 

1 .  Mode  terminology  has  been  clarified  and  a  new  framework  introduced.  The  “seeond  mode” 
is  not  really  a  seeond  mode  in  the  sense  of  a  discrete  spectrum  of  eigenvalues  and 
eigenfunctions  and  the  group  agreed  that  “fast”  mode  and  “slow”  mode  are  better  descriptions. 
However,  the  notion  of  seeond  mode  is  so  ingrained  in  the  community  it  may  be  hopeless  to 
try  to  affect  a  change.  A  meaningful  compromise  would  be  to  refer  to  the  seeond  mode  as  the 
Mack  mode  [7]. 

2.  A  tool  has  been  developed  for  the  projection  of  the  solution  of  the  linearized  Navier-Stokes 
equations  in  reacting  high  speed  boundary  layers  onto  diserete  modes. 

3 .  Employed  multi -mode  deeomposition  to  filter  out  the  unstable  mode  from  DNS  and  eompared 
the  result  with  a  theoretieal  prediction  based  on  the  method  of  multiple  seales  that  includes 
the  nonparallel  flow  effects  (in  collaboration  with  X.  Wang  and  X.  Zhong). 

4.  Utilized  PSE  approximation  to  study  the  discrete  mode  branching  (in  collaboration  with  Y. 
Eifshitz  and  D.  Degani  [25]  and  P.  Balakumar  [222]).  This  was  an  important  result  in  that  PSE 
will  always  stiek  with  the  correet  mode  and  the  diserete  branching  maybe  diffieult  to  interpret 
for  the  non-specialist. 

5.  Developed  boundary  layer  solvers  for  binary  mixtures  and  5-species  air. 

6.  A  code  for  stability  analysis  including  viscosity  and  real  gas  effeets  has  been  developed.  It 
was  shown  that  the  spectrum  with  real  gas  effects  is  similar  to  the  speetrum  in  a  calorically 
perfect  gas. 

7.  It  was  shown  that  an  inviscid  stability  analysis  captures  main  features  of  the  speetrum  and 
eigenfunctions.  The  real-gas  effects  considered  for  adiabatic,  non-catalytic  and  partially- 
catalytic  walls  could  be  interpreted  in  terms  of  the  wall-temperature  effect  for  the  seeond 
mode.  Behavior  of  the  mass  fraction  perturbation  is  governed  by  the  mean  flow  mass  fraetion 
gradient  and  the  vertical  velocity  perturbation.  Perturbation  of  the  chemical  source  term  is 
significant  for  the  second  mode  only  in  the  vicinity  of  the  wall,  and  it  results  in  a  change  in 
the  eigenvalues.  The  chemical  source  perturbation  plays  a  stabilizing  role,  and  the  meehanism 
is  governed  by  the  physics  in  the  vieinity  of  the  wall. 

8.  A  method  for  solving  receptivity  problems  in  presence  of  real-gas  effeets  has  been  formulated. 
The  considered  examples  indicate  that  there  is  an  effect  due  to  the  braneh  points  of  the  discrete 
speetrum  similar  to  that  in  the  calorieally  perfeet-gas  analysis. 
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9.  The  theoretical  model  of  the  second  mode  stabilization  by  porous  coatings  has  been  successfully 
validated  by  DNS  and  experiments  [3,  4,  8,  55,  61,  107].  Guidelines  for  the  choice  of  optimal 
porosity  parameters  have  been  formulated  [50].  The  porous  coating  laminar  flow  control 
technology  has  been  developed  to  the  readiness  level  suitable  for  flight  tests. 

10.  The  theoretical  model  [60]  of  boundary-layer  receptivity  to  thermal  fluctuations  has  been 
developed  to  estimate  for  the  upper  bound  of  transition  Reynolds  numbers  in  cases  of 
“absolutely  quiet”  free  stream  and  smooth  bodies,  where  the  forcing  source  is  associated  with 
indispensable  thermal  fluctuations  (Brownian  motion)  only. 

11.  Theoretical  models  of  supersonic  [9]  and  hypersonic  [56]  boundary-layer  receptivity  to  solid 
particulates  (such  as  dust  and  aerosols)  has  been  developed.  The  analytical  solutions  shed  light 
on  receptivity  mechanism  and  provide  foundation  for  prediction  of  particle-induced  transition 
in  high-altitude  flights,  where  particulates  can  be  the  major  forcing  source.  This  mechanism 
may  also  be  important  for  transition  in  wind  tunnels  and  ballistic  ranges  unless  special  filtering 
techniques  are  employed. 

12.  It  was  demonstrated  [62]  that  the  amplitude  method,  which  is  a  rational  physics-based  method 
for  transition  predictions,  can  be  implemented  for  certain  high-speed  configurations  including: 
T-S  or  second-mode  dominated  transition  induced  by  thermal  fluctuations;  second-mode 
dominated  transition  induced  by  solid  particulates;  C-F  dominated  transition  induced  by  small 
local  and/or  distributed  roughness  on  a  supersonic  swept  wing.  The  linear  receptivity  models 
were  developed  [9,  60,  63,  64]  to  evaluate  the  initial  spectrum  of  unstable  waves,  the 
asymptotic  theory  was  used  to  predict  the  downstream  propagation  of  boundary-layer 
disturbances,  and  empirical  amplitude  criteria  were  used  to  predict  the  transition  onset 
location.  It  was  shown  that,  in  the  foregoing  cases,  the  amplitude  method  can  be  realized  on 
economical  basis  with  keeping  the  numerical  requirements  nearly  the  same  as  for  the  e^ 
method.  This  effort  was  aimed  to  motivate  the  transition  study  groups  to  integrate  their 
knowledge  and  tools  in  a  rational  framework  of  the  amplitude  method. 

The  details  of  this  work  are  contained  in  the  Volume  VIII  supplement. 

3.4  Spatial  DNS  of  Breakdown  (Hermann  Fasel):  The  research  group  of  Fasel  investigated 
laminar-turbulent  transition  using  high-fidelity  direct  numerical  simulations.  The  objective  was  to 
conduct  DNS  in  close  collaboration  with  other  experimental  and  computational  efforts  to  identify 
relevant  mechanisms,  especially  the  nonlinear  breakdown  mechanisms  [54,  73] 

The  calculations  of  nonlinear  wave  interactions  [99,  153,  154,  155],  the  Gaster  wave  packet 
[149,150,  151,  152]  and  turbulent  spot  [148]  in  a  hypersonic  boundary  layer  were  remarkable 
because  of  the  similarities  with  incompressible  flow.  This  offers  some  promise  that  tackling  the 
breakdown  of  hypersonic  boundary  layers  is  not  as  formidable  as  once  thought.  It  makes  sense. 
The  second  mode  is  a  2-D  wave  as  is  a  T-S  wave.  The  nonlinear  interactions  with  background  or 
forced  3-D  modes  produces  a  lambda  vortex  structure  in  the  same  way  as  Klebanoff-type 
breakdown  in  incompressible  flow.  The  wall  streaks  have  been  observed  in  the  Purdue  Ludwig 
Tube  using  TSP  and  in  the  M6QT  at  A&M  using  IR  Thermography.  These  features  are  presently 
being  investigated  in  the  M6QT.  The  similarities  of  the  wave-packet  calculations  and  turbulent 
spot  calculations  with  the  incompressible  boundary  layer  are  very  encouraging.  It  seems  as  though 
we  have  a  foundational  starting  point  for  the  breakdown  of  hypersonic  boundary  layers. 
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For  modeling  the  porous  wall,  a  novel  immersed  boundary  method  was  developed  and 
implemented  into  the  in-house  developed  eompressible  Navier-Stokes  solver.  With  the  new 
scheme  it  was  possible  to  capture  the  effects  of  destabilization  through  very  narrow  porous 
coatings  as  predicted  by  theoretical  studies  by  Fedorov.  The  results  suggested  that  the  viscous 
dissipation  and  the  pressure  diffusion  play  a  major  role  in  the  stabilization/destabilization  process 
and  are  at  least  one  order  of  magnitude  larger  than  the  remaining  terms  computed  in  this  study.  In 
addition,  simulations  into  the  nonlinear  transition  regime  revealed  that  the  porous  wall 
considerably  delayed  the  breakdown  to  turbulence.  [49,73,74,75].  These  calculations  are 
complementary  to  Fedorov’s  work  in  section  (3.3). 

Linearized  Navier-Stokes  were  used  to  include  chemical  equilibrium  in  the  boundary-layer 
stability  calculations  [139]. 

The  details  of  this  work  are  contained  in  the  Volume  II  supplement. 

3.5  Surface  Roughness  and  Receptivity  (Xiaolin  Zhong):  In  addition  to  the  later  stages  of 
transition  work  of  Fasel,  two  other  areas  were  recognized  as  requiring  full  DNS  i.e.  surface 
roughness  and  receptivity.  These  are  summarized  in  the  Annual  Rev.  Fluid  Mech.  Article  [47] 

Developed  algorithms  for  surface  roughness  [5]. 

Examined  laminar  flow  control  using  2-D  roughness  [6,  10,  66,  67,  68] 
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